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Abstract 



We study non-adiabatic transitions in scattering theory for the time dependent 
molecular Schrodinger equation in the Born-Oppenheimer limit. We assume the elec- 
tron Hamiltonian has finitely many levels and consider the propagation of coherent 



X 



ctf 



states with high enough total energy. 

When two of the electronic levels are isolated from the rest of the electron Hamil- 
tonian's spectrum and display an avoided crossing, we compute the component of 
the nuclear wave function associated with the non-adiabatic transition that is gener- 
ated by propagation through the avoided crossing. This component is shown to be 
exponentially small in the square of the Born-Oppenheimer parameter, due to the 
Landau-Zener mechanism. It propagates asymptotically as a free Gaussian in the nu- 
clear variables, and its momentum is shifted. The total transition probability for this 
transition and the momentum shift are both larger than what one would expect from 
a naive approximation and energy conservation. 
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1 Introduction 



We study scattering theory for the time-dependent molecular Schrodinger equation 
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where the electronic hamiltonian h(x) is an m x m self-adjoint matrix that depends on the 
nuclear position variable x G M. The Born-Oppenheimer parameter e > denotes the fourth 
root of the electron mass divided by the mean nuclear mass. 

We compute the leading order asymptotics of nuclear wave functions associated with 
certain non-adiabatic transitions of the electrons. The Landau- Zener mechanism responsible 
for these makes them exponentially small in 1/e 2 as e — ► 0. 

Our most general result can be found in Theorem 15.11 Describing the most general 
situation requires the development of a significant amount of notation and some technical 
hypotheses. So, in this introduction, we describe two physically interesting special cases that 
illustrate the main consequences of our analysis in a simple situation. Theorems 16 . 1 1 and 16 . 21 
give precise statements of our results for these special cases. 

Suppose h(x) is a real 2x2 self-adjoint matrix that depends analytically on x and has 
limits /i(±oo) as x — > ±oo that are approached sufficiently rapidly. Denote the eigenvalues 
of h(x) by ej(x), and assume that e^ix) > ei(x) + 5 for all x G M, where 5 > 0. Near x = 0, 
assume e\ and e 2 have an avoided crossing, i.e., e2(x) — e\(x) ~ \/x 2 + S 2 close to x = 0, 
with 5 small but positive. Such an avoided crossing corresponds to complex crossing points 
zq and Zo, where the analytic continuations of e\ and e 2 satisfy ei(z ) = e 2 (zo), and z is 
close to the real axis, with z = 0(5). 

Let 0i (x) and ^(x) denote normalized, real eigenvectors associated with e\(x) and e2(x). 




Figure 1. A plot of typical electron energy levels involved in an avoided crossing. 



Among the nuclear wave functions we can accommodate are Gaussian coherent states 
that are defined by 



ipo(A, B, e 2 , a, i], x) 
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where the complex numbers A and B satisfy the normalization condition KeBA = 1. These 
states are localized in position near x — a, and in momentum near p = rj. Their position 
uncertainty is e\A\ and their momentum uncertainty is e\B\. For a thorough discussion of 
these wave packets, see |Hj. 

Choose E > snp x( -jj^ e2(x). For a state incoming from the left on the upper electonic level, 
choose r]_ > 0. We assume r\_ is large enough so that the classical energy 77^/2+e 2 (— oo) > E. 
There exists a solution to (fl.lj) whose large negative t asymptotics are given by 

where the nuclear part is a free Gaussian. Since the electronic levels are isolated from 
one another, the large positive t asymptotics of this solution are multiples of 02 (x), up 
to exponentially small errors in 1/e 2 . They have the leading behavior determined by the 
standard time-dependent Born-Oppenheimer approximation as e — > 0, see jHj: 

e i8 l{ e) e i(nt/2-e 2 (oo))t/e* ^^ + ^ ^ £ 2 ? ^ + ^ ^ ^ ^^ 

where e^ 1 ^ is some explicit phase, and the parameters Al, .Bi, Oi,77i are determined by the 
scattering properties of the classical Hamiltonian p 2 /2 + 62{x). 

Our interest lies with the leading order asymptotics of the non-adiabatic component of 
the wave function for large positive t and e — > 0. We prove in Theorem W>. II that these have 
the form 
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and we specify how the phase 9 + (e), the e-independent amplitude c > 0, the exponential 
decay rate a* > 0, and the parameters of the free Gaussian part A + , B + , a + , and r] + > 
are determined. As a corollary, the leading term of the transition amplitude A(e) (whose 
absolute square is the transition probability) is given by the quantity 

^( e ) = c e w+{€) e~ a ' /e \ as e -»• 0. (1.3) 

Let us describe the main features of this exponentially small transmitted part of the wave 
function. One may naively expect rj + to be determined by the energy conservation condition 

77 77 

— + e 2 (-oo) = y + ei(oo), 

but this yields the wrong value. The correct value is larger. Intuitively, this is due to the 
faster parts of the wave function behaving less adiabatically than the slower parts. Because 
this dependence on the speed appears in an exponent, it leads to an 0(1) change in the final 
momentum 77 + . In other words, the higher momentum components of the incoming state 
are much more likely to make a transition than the lower momentum components. Hence, 
after the transition, there are more fast pieces of the wave function, and the final average 
mometum is greater than one would naively expect from an energy conservation calulation 
based solely on the average incoming momentum. 

This also affects the transition amplitude which is larger than what is naively expected. 
It is asymptotically composed of an e-independent prefactor cq times an exponentially small 



quantity e~ a */ e , whose decay rate a* is related to that of the Landau-Zener decay rate for 
purely adiabatic problems. Actually, a* consists of the sum of the imaginary part of some 
action integral around the complex electronic eigenvalue crossing point zq and a contribu- 
tion that depends explicitly on the nuclear part of initial incoming state 1)1. 2j) . The action 



integral depends only on the electronic levels and reads /^ J2(E — e2(z)dz where ( is a loop 
in the complex plane based at the origin encircling zq. The contribution from the nuclear 
part of the wave packet depends on the shape of its momentum/energy density. It is that 
last contribution that makes the obvious candidate given by the imaginary part of the action 
integral taken at the classical energy E, miss the actual value of the decay rate a*. In that 
sense, (|1.3|) . which we could call a molecular Landau-Zener formula, cannot be determined 
from the usual adiabatic Landau-Zener formula with just the knowledge of the electronic 
levels and the classical nuclear momentum close to the avoided crossing. Indeed, our anal- 
ysis shows that we also need to take into account the details of the incoming wave packet 
to determine ([1.3)1 . This is why we resort to coherent states to get such accurate asymptotics. 

The way we obtain all our results is by employing a time-independent scattering theory 
approach that uses generalized eigenfunctions of the full Hamiltonian. We expand the wave 
function in terms of the generalized eigenfunctions and calculate the large \t\ asymptotics. 
For every incoming momentum k there is classical energy conservation, but a different proba- 
bility of making the non-adiabatic transition. We obtain the correct a* and r] + by computing 
the averages over k rather than by doing one calculation based on the average incoming mo- 
mentum ?7_. 

Remarks 

1. We obtain the analogous results when the incoming state is associated with the lower 
electronic level e x , provided that we keep the average total energy above both the levels. 

2. There are other components of the scattered wave function. For example, one should 
expect a reflected wave on the €2 electronic level and also a reflected wave on the e-y level. 
We prove that if the avoided crossing has a sufficiently small gap, then the other components 
are exponentially even smaller in 1/e 2 than the transmitted non-adiabatic term we compute. 

The second situation we describe in this introduction involves the same set-up as above, 
but with the Gaussian incoming states replaced by more general incoming coherent states. 
This example illustrates the second key feature that our analysis demonstrates: even if the 
incoming state is not Gaussian, the outgoing non-adiabatic transition state, generically, is 
Gaussian to leading order in e. 

For m = 1, 2, . . ., we define 

ip m {A, B, e 2 , a, r), x) = (1.4) 

2-/ 2 (m!) -i/2 ,4-™/* (A)^ Hm (^=£j MA, B, e 2 , a, V , x), 

where H m is the m th order Hermite polynomial. 
We now replace ([1.2)1 by 

e ^/2- e2 (-oo))iA 2 ^ A _ + lB _^ B 5 e 2 ) a _ + v ^ ^ x) ^ x y ( L5) 



Again, up to exponentially small errors, the large positive t asymptotics of the solution 
are multiples of 4>2{x). Their leading behavior is determined by the standard time-dependent 
Born-Oppenheimer approximation, 

e i6,{e) e i( V l/2-e 2 (oo))t/e* ^^ + -^ ^ ^ ^ + ^ ^ ^ ^^ 

where A±, B\, a±, r)i, and #i(e) are the same as in our first example. However, our Theorem 
16.21 shows that the leading order asymptotics of the non-adiabatic component of the wave 
function for large positive t again have the form of a freely propagating Gaussian 
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and display a pre-exponential factor of order e~ m . The values of a*, A + , B + , a + , and rj + are 
the same as in our first example, and we determine the pref actor c m . The numerics presented 
below clearly illustrate these features. 

Our most general result, Theorem 15 .![ extends these results in several ways. First, we can 
handle electron Hamiltonians h(x) that are m x m complex hermitian matrices which have 
two levels of interest that have an avoided crossing. These levels must stay well separated 
from the rest of the spectrum of h(x). Second, we can handle situations in which several 
levels display certain patterns of avoided crossings. For example, when two levels have an 
avoided crossing for one value of x, and one of those levels has another avoided crossing 
with a third level for some other value of x. However, in such cases, we can only study the 
non-adiabatic components for certain levels. The ones we can handle depend on the order in 
which the levels have the avoided crossings. Third, we can consider more general incoming 
states that do not have the form of the tpj's considered above. They are characterized by an 
energy (or momentum) distribution which is sharply peaked around some fixed energy, so 
that a semiclassical analysis can be performed. In such general cases also, the nuclear part 
of the non-adiabatic wave function is Gaussian and exponentially small, with a decay rate 
sharing the properties described above. 

The paper is organized as follows: In the rest of the Introduction, we review the relevant 
literature and present numerical results for the above examples. They show excellent agree- 
ment with our analysis. In Section 2, we set up the general problem we study. We state most 
of our hypotheses here and make precise the notion of avoided crossing. In Section 3, we 
study generalized eigenvectors of the full Hamiltonian. In particular, their WKB-type anal- 
ysis in the complex plane is performed here. We superimpose the generalized eigenvectors 
to generate solutions to the time-dependent Schrodinger equation and construct asymptotic 
scattering states in Section 4. Non-adiabatic transition asymptotics are studied in Section 
5, where our most general result is stated as Theorem 15 .11 Further properties and estimates 
on the energy and momentum shifts are provided in Section 5. Section 6 is devoted to the 
special case of interest where the nuclear part of the incoming state is a Gaussian or a Gaus- 
sian times a Hermite polynomial as in ()1.5|) . Finally, Section 7 contains the proofs of several 
technical results that are stated in the earlier sections. 

From this outline, one can see that our results depend crucially on the properties of 
generalized eigenvectors of the full Hamiltonian. We prove these properties by applying the 



ideas and results of Joye [Hj, [IS] that provide exponentially accurate WKB-type results in 
a generic avoided-crossing regime, generalizing earlier two-level adiabatic techniques from 
[T7j . [H], [IH]- See also [21], [2H] for stationary results of the same kind. That a complex 
WKB-type analysis plays an important role here should be no surprize. Indeed, in the ODE 
context of adiabatic-like problems dealt with in the references above, the complex WKB 
approach proved to be the most efficient method providing a quantitative analysis of the 
exponentially small leading order term of the Landau-Zener mechanism. See, however, J3] 
and [2] for a different successful approach of such problems, based on optimal truncation 
techniques. 

There are mathematical results on the exponentially small size of non-adiabatic tran- 
sitions in the Born-Oppenheimer approximation, and for related problems. See, e.g., [12J, 
|23j . [I], [21], [22]- However, to the best of our knowledge, there are no rigorous results on 
this topic in the literature that actually compute the leading asymptotics of non-adiabatic 
transitions in our time-dependent PDE setting. We have recently learned that Betz and 
Teufel, [3], are adapting techniques from [2] to the Born-Oppenheimer setup. They have for- 
mal and numerical results for specific electronic hamiltonians in agreement with ours. Also, 
rigorous results on the propagation of wave packets through avoided crossings, representing 
first attempts to unravel the molecular Landau-Zener mechanism, are obtained in ^H], [Hj- 
(See also (2EJ-) In those papers, the gap 5 shrinks to zero with e in such a way that the 
transitions are of order one, so that they can be computed by perturbation theory. This is in 
contrast to the present situation, in which S is small but fixed as e — > 0, and the transitions 
are exponentially small. 

Because of the importance of the Landau-Zener mechanism to molecular physics, there 
are relevant papers in the physics and chemistry literature. See, e.g., [I], [21], [2H]- 



1.1 Numerical Simulations for a Gaussian Initial State 

We now present graphical results of a numerical simulation in which the initial state is a 
Gaussian function associated with the upper energy level for a two level system. These plots 
are in very good agreement with the results of our analysis. 

We have numerically integrated equation (jl.lj) with e = 0.2 for the Hamiltonian function 



= 1 / 1 tanh(x) 

v ; 2 Vtanh(x) -1 



1 



The energy levels are ±- yl + tanh(x) 2 , and there is an avoided crossing at x = with 

a minimum gap of 1. The initial state is the eigenvector associated with the upper energy 
level times the Gaussian (fio(A + itB ,B ,e 2 ,r]t,ri,x), where A = B = r] = 1, with the 
initial time t = —10. The following two figures show the initial position and momentum 
probability densities, respectively. In both plots, the probablity of being on the lower energy 
level is zero. 




Figure 2. Position space plot at time t — — 10 of the probability density for being on the 
upper energy level (solid line), and 3 x 10 8 times the probability density for being on the 
lower energy level (dotted line). 




Figure 3. Momentum space plot at time t = —10 of the probability density for being on 
the upper energy level (solid line), and 3 x 10 8 times the probabilty density for being on the 
lower energy level (dotted line). 



The following two plots show the position and momentum probability densities at t — 9 
after the wave function has interacted with the avoided crossing. The component associated 
with the lower energy level has mean momentum 2.05. It is evident from the plot that it is 
greater than 2. 

The naive energy conservation calculation predicts the following: The total energy is 
E = rf/2 + 1/2 y / l + tanh(-10) 2 = 1.2071. After the transition to the lower surface, 

the kinetic energy should be this value plus n/2/2, so rJl/2 = 1.9142. This predicts a final 
momentum after the transition of r]i = 1.9566. 




Figure 4. Position space plot at time t = 9 of the probability density for being on the 
upper energy level (solid line), and 3 x 10 8 times the probabilty density for being on the 
lower energy level (dotted line). 




Figure 5. Momentum space plot at time t = 9 of the probability density for being on the 
upper energy level (solid line), and 3 x 10 8 times the probabilty density for being on the 
lower energy level (dotted line). 



1.2 Numerical Simulations for More General Initial States 

We next present the results for the same system as above, but where the initial Gaussian 
4>o has been replaced by 03. See (II. 4|) . Note that the transition amplitude is significantly 
larger than in the example above, and that the component of the wave function that makes 
the transition to the lower level is approximately a Gaussian. The value of epsilon e = 0.2 is 
not particularly small, so the component of the final state that does not make a transition is 
only approximately a 3 wave packet. We have chosen this relatively large value of epsilon 
to avoid numerical difficulties in integrating equation (jl.ll) . 

We should also note that the naive energy conservation calculation again predicts that 
the component of the wave function on the lower level should have mean momentum 1.9566. 
Since initial wave function has a greater momentum uncertainty than in the Gaussian exam- 
ple above, we see an even greater discrepancy between this prediction and the correct value. 
Our simulation yields a value of roughly 2.25. 




Figure 6. Position space plot at time t — — 10 of the probability density for being on the 
upper energy level (solid line), and 10 7 times the probabilty density for being on the lower 
energy level (dotted line). 
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Figure 7. Momentum space plot at time t — — 10 of the probability densityfor being on the 
upper energy level (solid line), and 10 7 times the probabilty density for being on the lower 
energy level (dotted line). 
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Figure 8. Position space plot at time t — 9 of the probability density for being on the upper 
energy level (solid line), and 10 7 times the probabilty density for being on the lower energy 
level (dotted line). 
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Figure 9. Momentum space plot at time t — 9 of the probability density for being on the 
upper energy level (solid line), and 10 7 times the probabilty density for being on the lower 
energy level (dotted line). 
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2 Hypotheses for the Electron Hamiltonian 

We begin with three general assumptions about the electron Hamiltonian h. We then impose 
two more assumptions that make precise the avoided crossing situations we can handle. 

HI: We assume z i— ► h(z) is a m x m matrix-valued analytic function that is analytic in 
z G p a = {z = x + iy : \y\ < a}, where a > 0. We assume h(z) is self-adjoint for z G M. 



Since we work in a scattering framework, we further assume: 
'2, c, and two matrices /i(±oo) 
sup \\h(x + iy) — h(±oo)\\ < 



H2: There exist v > 1/2, c, and two matrices /i(±oo), such that for all x G 1R, 

||//(.r+ it)) - //(±oc)|| < — ' 
\v\<<* 

where < x > denotes (1 + x 2 ) 1 ^ 2 . 



< x > 



2+v' 



The rate of convergence in this assumption can certainly be weakened. However, general 
scattering theory is not the main point of the present study. 

H3: We assume the spectrum a(h(x)) of h(x) consists of m non-degenerate eigenvalues 

a(h{xj) = {ej(x)} i= i,..., m , 

for any x G M U {±oo}. 



We let 



'iW> J 



1, • • • ,m, denote the corresponding eigenvectors, characterized up to 



constant phases by the following conditions 



ll&(*)ll = 1, 



and 



{(pj(x), <j>'Ax)) 



0, Vj 



•,m, 



(2-1) 
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-e,-(±oo)) = 


= 0{< x >- {2+u) ) 


and 






d n 
dx^^ x) 


-0i(±oo)) - 


= 0(<x>^ 2+ ^) 


for any n£N. 







where the the prime denotes the derivative with respect to x. The eigenvectors are analytic 
in some narrow open strip containing the real axis |20j . 

By using the Cauchy formula, it is easy to check that our hypotheses imply 

An 

(2.2) 



(2.3) 



We now make specific assumptions concerning avoided crossings for h. 

The idea is to assume h(x) belongs to a smooth family of electron Hamiltonians h(x, 5). 
When 5 = 0, we assume there are actual crossings. When 5 ^ 0, we assume there are no 
crossings for real values of x. The electron Hamiltonians we actually use have the form 
h(x, 5) for some small, but fixed value of 5. 

Our precise assumption is the following: 

H4: For each fixed 5 G [0, d], the matrix h(x,5) satisfies HI in a strip p a independent of 
5, and h(z,S) is C 2 as a function of the two variables (z,S) G p a x [0, d]. Moreover, h(-) 
satisfies H2 uniformly for 5 G [0, d], with limiting values h(±oo, 5) that are C 2 functions of 
5e[0,d\. 

Again, some of our results hold under weaker smoothness assumptions. 

We can deal with multiple avoided crossings, but cannot deal with all possible patterns 
of avoided crossings. The following assumption describes the ones we allow: 

H5: For each x G M and each 5 G [0, d], a(h(x, 5)) consists of m real eigenvalues 

a(h(x,5)) = {ei(x,(J), e 2 {x,8), ••-, e m {x,S)} c M. (2.4) 

When 5 > we assume these are distinct for x G [— oo, +oo] and are labeled by 

ei(x, 6) < e 2 (x, 5) < ■ ■ ■ < e m (x, 5). 

When 5 = 0, the eigenvalues are m analytic functions that have finitely many real crossings 
at x\ < X2 < ■ ■ ■ < x p , with p > 1. We assume the eigenvalues have m distinct limits as 
x — > — oo and as x — > oo. We label these eigenvalues ej(x, 0) in a way that is discontinuous 
in 5 near 5 = 0. This labeling is determined by the following conditions: 
i) For all x < x%, 

ei(x, 0) < e 2 {x, 0) < • ■ ■ < e m (x, 0). 

ii) For all j < I G {1, 2, • • • , n}, there exists at most one x r with 

ej(x r ,0) - e,(x r ,0) = 0, 
and if such an x r exists, we have 

— (e j (x r ,0)-e l (x r ,0)) > 0. (2.5) 
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iii) For all j G {1, 2, • • ■ , n}, the eigenvalue ej(x, 0) crosses eigenvalues whose indices are all 
superior to j or all inferior to j. 

Remarks: 

i) The parameter 5 can be understood as a coupling constant that controls the strength of 

the perturbation that lifts the degeneracies of h(x, 0) on the real axis. 

ii) Analyticity of the eigenvalues ej(x, 0) on the real axis follows from the self-adjointness of 

h(x,0). 

iii) The crossings are assumed to be generic in the sense that the derivatives of Cj — et are 

non-zero at the crossing x r . This ensures that when 5 > is small, the generic behavior 

()3.19|) holds at the corresponding complex crossing points. 

iv) When m = 2, H5 requires that the two eigenvalues have exactly one generic crossings 

when 5 = 0. 

v) The crossing points {x\, X2, ■ ■ ■ ,x p } need not be distinct, which is important when the 

Hamiltonian possesses symmetries. However, for each j = 1, • ■ ■ ,n, the eigenvalue ej(x,S) 

experiences avoided crossings with e_j + i(x, 5) and/or ej_i(x, 5) at a subset of distinct points 

\X T1 , • ■ • , X r - j <— \X\, X2, ■ , Xpj. 

For certain results, we also impose the condition that these avoided crossings be generic 
in the sense of [7j and (T3J. This condition essentially says that the low order Taylor series 
coefficients of certain quantities do not vanish at the crossing when 5 = 0. 

H6: Near an avoided crossing of ej(x, 5) and e n (x, 5), there exist a > 0, b > 0, and c G M, 
such that 



e n (x,5) - ej(x,6) = ± ^ax 2 + 2cx5 + b 2 5 2 + R 3 {x,5), (2.6) 

where c 2 < a 2 b 2 and #3(2;, 5) is a remainder of order 3 in (x, 5) close to (0, 0). 

Our final hypothesis involves both the electron Hamiltonian and an interval of energies, 
A. We ultimately consider states of the full Hamiltonian whose energy is concentrated in 
A, with A high enough that scattering onto all the electron energy levels is possible. An 
energy range that satisfies this condition can always be chosen for some strip p a , provided 
the minimum value in A is large enough. 

H7: The interval A G M is compact and has non-empty interior. Furthermore, it is chosen 
so that 

inf \E-eAz,8)\ > 0. 

EGA J V 

5S[0,<5] 

3 Generalized Eigenvectors 

For energies E G A, we construct generalized eigenvectors for the full Hamiltonian. For the 
time being, the parameter 5 > is fixed and we drop it in the notation. The generalized 
eigenvectors are solutions ^(x,E,e) G (D m to the time-independent Schrodinger equation 

( ~ 1 & 2 + Hx) ) * (X ' E} e) = E * (X ' E ' e) - (3 ' 1} 

12 



For each E G A, the set of such solutions is 2m-dimensional, and individual solutions can 
be characterized by their asymptotics at x = — oo (or at x = oo). 
Let 



$(x,E,e) = I . „ „\;., : J6 (T 2m . 



ie 2 f ^(^^e) 



*e 2 |-$(x,£,e) = H(x,E)${x,E,e), (3.2) 



Then (J3.1J) is equivalent to 
<9x 



where 

/ Q E \ 

H(x,E) = € M 2m ((D) and 

EE-h(x) > 0, for all x€iRU{±oo}. (3.3) 

Here, E denotes the identity matrix in (T m . Note that the matrix H (x, E) is not self-adjoint, 
but satisfies the relation 

H(x,E) = J H*{x,E) J, where J = (% jT\ . (3.4) 

The small e asymptotics of solutions to ()3.2|) are studied in [T3] . Of particular importance 
to us is Section 7 of JT^], which is devoted to the computation of exponentially small elements 
of the related ^-matrix that we describe below. We apply the results of [TH] to (J3.2J) . keeping 
track of the dependence on E. 

By our hypotheses on E and h(x), the spectrum of H(E, x) consists of 1m distinct real 
eigenvalues 

a{H{x,E)) = {k^x,E)} T =+:~ m , with 



k](x,E) = rkj(x,E) = T^2(E-e j (x)) e JR. (3.5) 

Note that the fcj's correspond to the classical momenta associated with the classical poten- 
tials ej(x). 

A set of corresponding eigenvectors { Xj{ x > E) } * s given (in block notation) by 

**<*■*> - { k *(XJ e ^ (3 - 6) 

From these we produce new eigenvectors 

^(x,E) = X ]{x,E) \ (3.7) 

that satisfy the normalization convention (|3.1(J|) below, that wasadopted in [T5\. This nor- 
malization is motivated by the following: We can write 

H(x,E) = £ k^(x,E)Pj(x,E), 
13 



where {PUx, E)} denotes a set of non-orthogonal projections onto the eigenspaces of H (x, E). 
If we define 

1 / 7 (x) \ 

e ' (x ' E) -A^L*))> (3 - 8) 

then it is easy to check that 

PJ(x,E) = \ X ]{x,E)){e]{x,E)\, (3.9) 

where we have used the bra-ket notation relative to the scalar product in <D 2m . We use the 
same notation for scalar products in (D m and (D 2m , since no confusion should arise. 
We now see that the eigenvectors (|3.7|) satisfy the normalization conditions 

PJ(x,E) £^{x,E) = 0, and 

(3.10) 
(^J(0,£), J^(0,E)) = r G {-1,1}. 

We note that H, fcj, xj, PJ, and <£>J are analytic functions of x and E when these 
variables are in a neighborhood of M x A. More precisely, if A = [.Ei,^], we define 
Dp = {z G (U : dist(z, A) < /3}, and these functions are analytic in p a x Dp, for a and (3 
small enough. Here a must be chosen small enough so that Cj and <pj are analytic in p a , (see 
|2Ti]). and /3 must be small enough so that \E — ej(x)\ > in p a x Z}^. 

We later make use of larger values of a in order to take advantage of the generic multi- 
valuedness of e, and 4>j as functions of x. 

From J3], we now see that any solution to (J3.2J1 can be written as 

$(x,£,e) = Y, $(x,E,e) e-' X W^*"^ <pT(x,E), (3.11) 

where the scalar coefficients d" G (F satisfy the equation 



<9x 
with 



£- cj(x, £, e) = 53 aj^x,^) e'^^^)-^^))^ 2 c -(x,E,e), (3.12) 

a jl \ x i^) 



J || (pj{x,E) |p 

We can rewrite ()3.12)1 as an integral equation 

c T j(x,E,e) = cJ(ar ,J5,e) 



+ / X £ a^(x',E)e i S^' rk ^^- ffk ^'^ d y^ c?(x',E,e)dx'. (3.13) 



' x ° J,ff 
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As we shall soon see, our hypotheses imply the existence of the limits lim x ±oo cj(x, E, e) 
cJ(±oo, E, e), so that with the notation 



{ c[{x,E,e)\ 



e <D r 



c T (x,E,e) = 

we can define an associated S'-matrix, S G M 2m ((Z7), by the identity 
S(Ee) ( c+ (-°°> E > € )\ (c + (+oo,E,e)\ 

This S'-matrix naturally takes the block form 

\S-+(E,e) S—{E,e))' 

Due to the symmetry (|3.4|) . it also satisfies the relation (see [15J), 



S(E,e) 



(3.14) 



(3.15) 



S-\E,e) = RS*{E,e)R, where 



R 



E Q 
(D -E 



Its elements describe transmission and reflection coefficients at fixed energy E which play key 
roles in our analysis. The off-diagonal elements are exponentially small and their asymptotics 
are determined in [To] . 

With this notation, the generalized eigenvectors are given by 



*(x,E,e) 



E 



2kj(x,E) 



<j>j{x) 



(3.16) 



x 



cUx,E,e)e- l fo k ^ E)dy/e2 + cj(x, E, e) e l So k ' 



(V,E) dy/e 



Since Jq kj(y,E)dy ~ x kj(±oo, E) = x J2 [E — ej(±oo)) as x — > ±cxd, the component 
of (|3.16|) that describes a wave traveling from the left to the right is labeled by — , and the 
component that describes a wave traveling from the right to the left is labeled by +. Note 
also that ()3.16|) is simply a WKB decomposition of the generalized eigenvectors. 

We now state some of the general properties of the coefficients dj(x,E,e) and of the 
phases e l Jo j ( y ' > y ' e that allow us to justify the scattering results described above. 

Lemma 3.1 Our hypotheses on h(x) imply the following, uniformly for E G A and all 
n G N: 



< C x (n) < 



gr. 



dE> 



kj(x,E) < 62(71) < 00, 



0" 



BE 



- (k](x, E) - fcJ(±oo, E)) = 0(< x > 

15 



-(2+«/)> 



and (3-17) 

as x — > ±00. (3.18) 



Thus, if we define co>J(±oo, E) = / °° (kj(y, E) — &;J(±oo, E)) dy, we further have 
f k](y,E)dy = xk°(±oo,E) + uJ(±oo,E) + r?(±,x,E) 

•/ U 

where, uniformly in E and for all n G N, 

' )'lt 

— r?(±,x,E) = 0(<x>-^ 1+ % as x - ±00. 

Moreover, the limits cJ(±oo, E 1 , e) asi-> ±00 exisi, and as |x| — > oo ; 

<9" 



<(*,£, e) = 0(1), for n = 0, 1, 



uniformly for E G A. ^/so ; as 1^ ±00 and uniformly for E G A, 

cj(x,£,e) - cJ(±oo,£,e) = 0(< s >" (1+I/) ), and 

<9 



5E 



c°(x,E } e) - dj (±00, E, e)\ = 0(< x > 



Remarks: 

1. This lemma is proved in Section [7| 

2. Some error terms in the lemma may depend on e in a singular way as e —>■ 0. This will 
not matter because they will vanish in the scattering framework we adopt. 



3.1 Complex WKB Analysis 

All the information about transmissions and transitions among the asymptotic eigenstates of 
the electronic Hamiltonian is contained in the asymptotic values of the coefficients cj(x, E, ±00) 
defined by (|3.13|) . and hence, in the matrix S(E, e). We extract this information by mimick- 
ing the complex WKB method of J3] , while keeping track of the E dependence. 

The complex WKB method requires hypotheses on the behavior in the complex plane of 
the so-called Stokes lines for the equation (J3.2J1 in order to provide the required asymptotics. 
These hypotheses are global in nature, and in general, are extremely difficult to check. 
However, in the physically interesting situation of "avoided crossings," they can be easily 
checked. We restrict our attention to these avoided crossing situations that are described 
below. 

We consider the coefficients Cj that are uniquely defined by the conditions 

cj(— 00, E, e) = 1, and c£(— 00, E, e) = 0, for all (k, a) 7^ (j, r). 

The key to the complex WKB method lies in the multivaluedness of the eigenvalues and 
eigenvectors of the analytic generator H(x, E) of (|3.2|) in the complex x plane. For any fixed 
E G A, H (-, E) is analytic in p a , and the solutions (|3.11|) to (J3.2J1 are analytic in x as well. 
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However, the eigenvalues and eigenvectors may have branch points in p a whose properties 
are inherited from those of the eigenvalues and eigenvectors of h(-). 

Analytic perturbation theory as described in [23] states that the eigenvalues and eigen- 
projections of h(x) for real x are analytic on the real axis and admit analytic multivalued 
extensions to p a . The analytic continuations of the eigenvalues have branch points that are 
located on a set of crossing points 

fi = {z G p a \M : ej(z ) = efc(z ) for some j, k and some analytic continuations}. 

Recall that for 5 = 0, the eigenvalues are analytic at any crossing points on the real axis. 
This follows from the self-adjointness of h(-) on the real axis. Note also that Q = Q by the 
Schwarz reflection principle. 

Generically, at a complex crossing point Zq G Q, we have the following local behavior, 
where c G (D is some constant 

ej (z) - e k (z) ~ c(z-zo) 1 / 2 (1 + O(z-Zo)). (3.19) 

The eigenprojectors of h(x) also admit multivalued extensions in p a \ Q, but they diverge at 
generic eigenvalue crossing points. We only have to deal with generic crossing points. 

To see what happens to a multivalued function / in p a \Q when we turn around a crossing 
point, we adopt the following convention: We denote by f(z) the analytic continuation of / 
defined in a neighborhood of the origin along some path from to z. Then we perform the 
analytic continuation of f(z) along a negatively oriented loop that surrounds only one point 
z G Q. We denote by f(z) the function we get by coming back to the original point z. 

We define Co to be a negatively oriented loop, based at the origin, that encircles only z 
when Im z > 0. When Im 2 < 0, we choose Co to be positively oriented. 

We now fix zq with Im z > 0. If we analytically continue the set of eigenvalues \e,j (z) }" 
along a negatively oriented loop around z G Q, we get the set {ej(z)}" l =1 with 



m 
1; 



ej(z) = e noU) (z), for j = 1, • • • , m, 

where 

7r : {l,2,...,m}->{l,2,--.,m} (3.20) 

is a permutation that depends on z . As a consequence, the eigenvectors (|2.1|) possess 
multivalued analytic extensions in p a \Q. The analytic continuation (f)j(z) of <j>j{z) along 
a negatively oriented loop around zq G f2, must be proportional to <frir (j)(z). Thus, for 
j = 1, 2, • ■ ■ , m, there exists #j(Co) £ ®, such that 

h{z) = e- i6 ^^ Q{j) {z). (3.21) 

We now turn from h(x) to H(x, E). 

From Hypothesis H7, ()3.5|) . and ()3.7|) . we see that the set of crossing points for the 
eigenvalues ±kj(x, E) of H(x, E) is independent of E and coincides with Q. 
Moreover, for j = 1, • • • , m, we have 

k]{z,E) = K 0{j) (z,E), (p]{z,E) = e -^°)^ o(i) (*,£), 
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where the prefactor e~ l9: >^ is independent of E. 

The above implies a key identity for the analytic extensions of the coefficients cj(z, E, e), 
z G p a \£l. Since the solutions to ()3.2|) are analytic for all z G p a , the coefficients cj must 
also be multi- valued. In our setting, Lemma 3.1 of jTHJ implies the following lemma. 

Lemma 3.2 For any j = 1, • • • , m and r = +, — , we have 

^(z,E,e) e ~ i k k >' E)du/e2 e~ ie ^o) = ^(^ £, e ) (3.22) 

where (o, $j(Co) an d 7To(i) ar & defined as above and are independent of E G A. 



Remark: Since Q is finite, it is straightforward to generalize the study of the analytic 
continuations around one crossing point to analytic continuations around several crossing 
points. The loop (o can be rewritten as a concatenation of finitely many individual loops, each 
encircling only one point of Q. The permutation ttq is given by the composition of associated 
permutations. The factors e~ tdj ^ ' in (J3.21)) are given by the product of the factors associated 
with the individual loops. The same is true for the factors exp ( — i f^ Q kj(z, E) dz/e 2 J in 
Lemma 13.21 

We now describe how we use the above properties. The details may be found in [TH] . 

The idea is to integrate the integral equation corresponding to (J3.13)) along paths that 
go above (or below) one or several crossing points, and then to compare the result with the 
integration performed along the real axis. As z — > — oo in p a these paths become parallel 
to the real axis so that the coefficients take the same asymptotic value c T m {— oo, E, e) along 
the real axis and the integration paths. Since the solutions to (J3.2J) are analytic, the results 
of these integrations must agree as Re z — > oo. Therefore, (J3.22J) taken at z = oo yields 
the asymptotics of ^w-Joo, E, e), provided we can control Cj(z,E,e) in the complex plane. 
We argue below that this can be done in the so-called dissipative domains (See 0, |S]), as 
proven in [T3]. We do not go into the details of these notions because another result of [T3] 
will enable us to get sufficicient control on Cj(z,E,e) in the avoided crossing situation, to 
which we restrict our attention. 

We define 

A™(x,E) = J o (k'(y,E) - k](y,E)) dy. 

By explicit computation, using formula ()7.3|) in (j3.13|) . we check that (J3.13J) can be extended 
to p a \ Q. We next integrate by parts in (|3.13|) . to see that (|3.13|) with x = — oo can be 
rewritten as 



c^(z,£,e) = «5, m 5 T _-*e 2 y> " mZ ^f ; j™*™* %(z, E, e) (3.23) 

J b k^ m (z,E)-kf(z,E) l{ } K } 



a ml\ Z ->E) iA T,7 ,(z,E)/e 

— g ITll \ ' '< 

,E)-kf{z,E) 

+ ie 2 V f (— - ^!^l3 ^ e ^(^)A 2 zUz' E e) dz' 

t- J-oo \ dz' %T. (z>. E) - k?(z' E) K ' ' ; " 



r* aZ(z>, E) a% (z>, E) ^ p{AE)/e , ^ , , 

,^J-ockT(z'E)-k?(z>E) pV ; 



,^ >9 J-< x >kJ n (z',E)-kf(z' : E) 



as long as the chosen path of integration does not meet Q. Here, ~ denotes the analytic 
continuation along the chosen path of integration of the corresponding function defined 
originally on the real axis. This distinguishes c^(oo, E, e) from c^(oo, E, e) computed along 
the real axis as x — > oo. These quantities may differ since the integration path may pass 
above (or below) points of Q. 

If the exponentials in ()3.23|) are all uniformly bounded, as it is the case when the inte- 
gration path coincides with the real axis, it is straighforward to get bounds of the type 

5 T m (z,E,e) = 5 jm 5 T . + E (e 2 ). (3.24) 

In our context, all quantities depend on £ e A. However, by mimicking the proof of 
Proposition 4.1 of 15], it is not difficult to check that the estimate (J3.24)) is uniform for 
E G A. For later purposes, we note that by differentiating (|3.23|) . f^c^z, E, e) is uniformly 
bounded for < e < eo and E G A for any fixed €q. 

Again, as is well known, the existence of paths from — oo to +oo along which the expo- 
nentials do not blow up and which pass above (or below) points in Q is difficult to check in 
general. It is linked to the global behavior of the Stokes lines of the problem. See e.g., |Sj, 
[5]. This property goes under the name "existence of dissipative domains" in |15j . 

We avoid these complications by restricting attention to avoided crossing situations where 
the existence has been proven [T5] . 

When dissipative domains exist, ()3.22|) and ()3.24|) imply 



-"iro(j) 



(oc,£,e) = e ~ l ^ k ^ E)du,e2 e~ ie ^ (1 + O^e 2 )), (3.25) 



where the Os(e 2 ) estimate is uniform for E G A. This is the main result of Proposition 4.1 
in [13] in our context, under the assumption that dissipative domains exist. 



3.2 Avoided Crossings 

We now explore the avoided crossing situation, alluded to above, that allows us to avoid 
considerations of the dissipative domains. We now assume that h(x) has the form h(x, 5) 
and satisfies Hypotheses H4 and H5. 

We first check that the allowed pattern of avoided crossings for a(h(x, 5)) can be trans- 
fered to the eigenvalues of H(x, E, 5), obtained from h(x, 5) by (|3.3|) . 

From the explicit formulae (|3.5|) , we see immediately that x c G M is a real crossing point 
for the eigenvalues e,(x, 0) and e/(x, 0) of h(x, 0) if and only if it is a real crossing point for 
the analytic eigenvalues kj(x, E, 0) and kf(x, E, 0) of H(x, E, 0), for r = +, — . Moreover, 



£(e l (x,0)-e J (x,0)) 



^-(k^(x,E)-kl(x,E)) 

V ^ — ^ c 

so that the real crossings for H(x, E, 0) are also generic, in the sense of ()2.5|) . 
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Remark: Our assumption H7 on the parameter E forbids real crossings between eigenval- 
ues k~j(x, E, 0) and kf (x, E, 0), with a ^ r. 

Regarding the ordering of the eigenvalues of H(x, E, 5), if those of h(x, S) are ordered as 
in (J2.4J) . we have 

- k t (x, E,S) <■■■ < -k m (x, E, 5) < < k m (x, E,S) <■■■ < k x {x, E, 5). (3.26) 

This means that the pattern of the crossings for the group of eigenvalues {— kj(x, E, 0)}j = i... jm 
is the same as that for the eigenvalues {ej(x, 0)} J= i 5 ... ?m . The pattern of the crossings for the 
group {kj(x, E, 0)}j = i t ... tTn is the reflection with respect to the horizontal axis of the one for 
{ej(x, 0)}j=i l ... tm . Therefore, assumptions H5, i), ii), hi) are also satisfied for the eigenvalues 
of H (x, E, 0), for a relabeling from 1 to 2m of ()3.26|) with 5 = 0, and x close to — oo. 

To any given pattern of real crossings for the eigenvalues {ej(x, 0)}j=i,..., m of h(x, 0), we 
associate a permutation 7r of {1,2, • • • , m} as follows. Assume the eigenvalues are labeled 
in ascending order at x = — oo, as in property i) of H5. If e.,(oo, 0) is the k th eigenvalue in 
ascending order at x = oo, the permutation tt is defined by 

tt(j) = k. (3.27) 

We call ir the permutation associated with a(h(x, 0)). For small 5 > 0, the real crossings 
turn into avoided crossings on the real axis and conjugate complex crossing points appear 
close to the real axis. Then tt corresponds to the permutation ttq (J3.20J) for a loop Co that 
surrounds all complex crossing points in the upper half plane that are associated with the 
avoided crossings. 

These properties of corresponding patterns of real crossings of the spectra of h(x, S) and 
H(x, E, 0) immediately yield the following convenient relation between the permutation n 
associated with a(h(x, 0)) and the permutation II associated with a(H(x, E, 0)). If we denote 
II by the obvious notation 

U(j,r) = (k,a), 

then we have 

n(j,r) = (7r(j),r), for all (j, r) G {1, ■ • • , m} x {-, +}. 



We can now restate the main result of [12] that describes the asymptotics of the coef- 
ficients defined in (J3.13)) . adapted to our scattering framework for incoming states entering 
from the left. (See (|3.1fjjl .) Intuitively, this result says that for small 6 > 0, dissipative do- 
mains exist, provided the pattern of real crossings satisfies H5. Therefore, estimates of the 
type ()3.25|) are true for certain indices j and n, determined by the permutation (J3.27J) . It is 
not difficult to see that the permutation ir describes the successive exchanges of eigenvalues 
one gets by following a path in the complex plane that goes above or below all complex 
crossing points of the eigenvalues Cj (x, 5) that are associated with the avoided crossings. 
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Theorem 3.1 Let h(x, 5) satisfy H4 and H5. If 5 > is small enough, the 7t(j),j elements 
of the matrix S (E, e) in \3.15\i . with n(j) defined in \3.27§ , have small e asymptotics for 
all j = 1, • • • , m, given by 

S- m {E, e) = n e- im > 5) ^ hW)d ^ (l + E Ae 2 )) , <j) > \ 

1=3 y < j 

where, for n(j) > j (respectively n(j) < j), b,l = j,---, rc(j) - 1 (resp. I = j, ■ ■ ■ , n(j) + I), 
denotes a negatively (resp. positively) oriented loop based at the origin which encircles the 
complex crossing point z r only (resp. z^) corresponding to the avoided crossing between 
ei(x, 5) and ei + \(x, 5) (resp. ej_i(a;, S)) at x r . The L ki(z, E, 5) dz denotes the integral along 
(i of the analytic continuation ofki(0,E,S), and 9i(Ci) is the corresponding factor defined by 

(EH. 

Remarks: 

i) Revisiting the proof of this theorem in [TS|, we see that we can choose 5 > small enough 
so that dissipative domains can be constructed uniformly for E G A. This stems from the 
formula 



kj(x,E,0)-ki(x,E,0) 



2(e l (x,0)-e j (x,0)) 



k j (x,E,0) + ki{x,E,Oy 

whose denominator can be controlled, close to the real axis, uniformly for E G A. 

ii) When there is only one avoided crossing between level j and j + 1 stemming from a real 

crossing at x = xq, we have j + 1 = ir(j) . The theorem says 

S- u - +1)>J (E,e) = e-^e^ mdz " (l + E ^)) , 

where the negatively oriented loop Q encircles the corresponding complex crossing point 
Zq, with Im zq > 0. Similarly, interchanging the roles of j and j + 1, it yields with Q the 
conjugate of the loop Q, 

S Iu+l) (E,e) = e-«.^e i ^ t ''- E " d "'\l + E . t ^)). 
iii) Since the eigenvalues are continuous at the complex crossings, we have 

lim Im / kj(z, E, 6) dz = 0, for all j — 1, • • • , p. 

It is shown in [T3] that 

lim Im 9j(Q, 5) = for all j = 1, • • • ,p. 

iv) The Oe,s(^) errors in Theorem 13 . II depend on 5, but it should be possible to get estimates 
which are valid as both e and S tend to zero, in the spirit of ^3], J2I], and [25] . 
v) This result shows that at least one off-diagonal element per column of the S-matrix can 
be computed asymptotically. However, it is often possible to get more elements by making 
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use of symmetries of the S- matrix. See ^3] and ^H] • 



In our avoided crossings context, transitions of the coefficients between states that cor- 
respond to electronic levels that do not display avoided crossings, i.e., that are separated 
by a gap of order 1 as 8 — ► 0, are expected to be exponentially smaller than the transitions 
we control by means of Theorem 13.11 as 8 shrinks to zero. Since the coefficients in the ex- 
ponential decay rates given by the theorem vanish in the limit 8 — > 0, it is enough to show 
that the decay rates of the exponentially small transitions between well separated levels are 
independent of 8. 

That is the meaning of the following proposition, which draws heavily upon JH] and J3] 
and is proven in Section [7| 

Proposition 3.1 Let F(x,5) be a n x n matrix that satisfies H4 ; except for the condition 
that F(-,8) be self- adjoint. Suppose its eigenvalues {fj(x,$)}j=i,-..,m that satisfy H5. Further 
assume that the eigenvalues can be separated into two groups &\{x, 8) and o~2{x, 8) that display 
no avoided crossing, i.e., such that 

inf dist(o"i(x, 8), cr 2 (x, 8)) > g > 0. 

x£paU{±oc} 

Let P(x, 8) and Q(x, 8) = E— P(x, 8) be the projectors onto the spectral subspaces corre- 
sponding to o"i(x, 8) and o~2{x,8) respectively, and let U e (x,xo,8) be the evolution operator 
corresponding to the equation 

ie 2 — UJx, xq, 8) = F(x,5) UJx,xq,5), with UJxq, xo, S) —E (3.28) 

ax 

Then, for any 5 > 0, there exists €q(5), C(S) > depending on 5, and T > independent of 
5, such that for all e < e (S), 

lim \\P(x,5)U e (x,x ,5)Q(x ,5)\\ < C(5)e- r ^. 



x — >oo 

XQ — > — OO 



Remark: This proposition implies that reflections, i.e., the transitions from wave pack- 
ets traveling to the right to wave packets traveling to the left, on any electronic level, are 
exponentially smaller than transitions associated with the avoided crossings in which the 
propagation direction is not changed. This is a consequence of Hypothesis H7 which implies 
that complex crossings between fc+ and kf , are far from the real axis for any j,l 6 {1, • • ■ , m}. 

Let us investigate more closely the analytic structure of kj(z, E, 8) in our avoided crossing 
regime characterized by H4 and H5, in order to deduce the properties of the exponential 
decay rates Im J^. kj(z,E,8) dz. We do so for the k/s that correspond to electronic eigen- 
values ej(z, 8) and e n (z, 8) that experience only one avoided crossing, i.e., we take n = j ±1. 
We can thus drop the index j in Q in the notation. We follow [T^j where a similar analysis 
is performed, sometimes refering to results proven there. The general case is dealt with by 
adding the corresponding contributions stemming from each individual avoided crossing. 
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We can assume that the avoided crossing takes place at x — 0, i.e., 

e,(0,0) = e„,(0,0) = e c , 

where e c is the electronic eigenvalue at the crossing. We also define the momentum k c (E) at 
the crossing point by 



k c (E) = p(E-e c ) 
and the quantity T Q (S) by 

r (6) = Im I e j (z,5)dz . (3.29) 

This quantity is the exponential decay rate given by the Landau-Zener Formula for a (time 
dependent) adiabatic problem with hamiltonian h(t,S). See jHj. In Section [7| we prove 

Lemma 3.3 With the above notation, we have the following as 5 — > 0, uniformly for E G A ; 
Im jf y/2(E- ej (z,S))dz = ^ + 0(6"), 
^ Im jf y/2(E- ej (z,5))dz = - ^| + 0(5 3 ), and 

■^2 Im l y/2(E-e s (z,6))dz = 3 ^ + 0(5 3 ), 

where < r (<J) = 0(5 2 ). 

Remarks: 

i) This implies that Im L J2(E — 6j(z, 5)) dz is a positive, decreasing, convex function of 
E on A. This remains true when the transition is mediated by several avoided crossings, 
ii) The first relation can be interpreted as saying that in our Born-Oppenheimer context, the 
(time dependent adiabatic) Landau-Zener decay rate at fixed energy E has to be modified 
in order to take into account the classical velocity k c (E) at the crossing, 
iii) More precise estimates will be derived below, further assuming H6. 



4 Exact Solutions to the Time— Dependent Schrodinger 
Equation 

We now construct solutions to (jl.lj) by taking time-dependent superpositions of the gen- 
eralized eigenvectors ^(x,E,e), where E G A. These superpositions depend on an energy 
density Q(E,e) that can be complex and may or may not depend on e. We always assume 
that the following condition holds: 

CO : The density Q(E, e) is C l as a function of E G A, for fixed e > 0. 
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In this Section, the parameter 5 > is kept fixed and plays no role. We therefore drop 
it from the notation and work under hypotheses HI, H2, and H3. 

We define 

ij>(x,t,e) = f Q(E,e)ty(x,E,e)e- ltE/e2 dE 

J A 

£ 0j(x) I W E > e) c°(x, E, e) e-'i: *^)'W' a e ~ UE ^ dE 
j=i,..., m , <r=± Ja J2kj(x, E) 



J2 i%(x,t,e). (4.1) 

j=l,---,m <r=± 

Here ipj asymptotically describes the piece of the solution that lives on the electronic state 
4>j and propagates in the direction characterized by a. Since the integrand is smooth and A 
is compact, ip(x,t,e) is an exact solution to the time-dependent Schrodinger equation (jl.lj) . 
Note that this solution is not necessarily normalized. 

The following lemma, whose proof can be found in Section Q gives a bound that we use 
to understand the large t behavior of i/jJ(x, t, e). It is a simple corollary that the state (|4.1jl 
belongs to L 2 (M). 

Lemma 4.1 Assume HI, H2, H3 and CO. Let 



and 



K + = sup kj(aoo,E) < oo 

j=l,—,m EeA, a=± 



K_ = inf kAaoo.E) > 0. 

j=l,-,m EGA, o=± J 



Fix a G (0, 1). Then, for either t = 0, or for any x ^ and t ^ 0, such that 

\x/t\ > K+/(l-a) or \x/t\ < !T_/(l + a), 

we have 

ipj(x,t,e) < C e /\x\, with C e independent of t, 

where the estimate is in the norm on(E m . 

We now introduce freely propagating states ip(x, t, e, ±) e L 2 (M, (C m ) that describe the 
asymptotics of the solutions ip(x, t, e) as t — > ±oo. We use these asymptotic states when we 
study the scattering matrix for (jl.lj) . We let 

i/>(x,t,e,±) (4.2) 

t Q(E,e) e _ UE/e 2 E ^ _ tixkJ(± ^ E)+ ^ {± ^ E) y e , dE 



X ^ x 



i=l,...,m «r=± ' JA ^2kj(±00,E) 

X i>J(x,t,e,±) 
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These states are linear combinations of products of free scalar wave packets in constant 
scalar potentials times eigenvectors of the electronic Hamiltonian. Their propagation is thus 
governed by the various channel Hamiltonians. 

Proposition 4.1 Assume HI, H2 ; H3 and CO. In L 2 (M) norm as t — > ±oo ; we have 

U(x,t,e) -4(x,t,e,±)\\ = O e (l/\t\). (4.3) 

Remarks: 

i) The estimate (J4.3J1 depends on e. 

ii) By a change of variables, we immediately obtain the following corollary. 

Corollary 4.1 The density of the component of the asymptotic momentum space wave func- 
tion on the j th electronic level as t — > ±oo is 



Ik 



a J-Q(E(k),e) c°(±oc,E(k),e) e 



-iui?(±oo,E(k))/e 



2 



Here E{k) = k 2 /2 + ej(±oo) and a = —/+ for waves traveling in the positive/negative 
direction, respectively. 

iii) Consider a solution i[)(x,t,e) traveling in the positive direction and associated with the 
electronic eigenstate <pj in the remote past. It is characterized by c%(— oo, E, e) = 5kjo~ a -, 
and as t — > — oo, it is asymptotic to 

#M,e,-) = Mx) I Q ^ E > 6 ) e -^A 2 ^(^(-^^-(-oo,^)/^ djg , (4.4) 



As t — > +oo, the component of this state that has made the transition from state j to state 
n is asymptotic to the vector ip~(x, t, e, +). It is given in terms of the matrix S by 

n ( x ) / gffi e ) e -^A 2 5-r(E, e) e +«(*M+°°.£)-"n (+°°,£))A 2 rfE- (4 _ 5) 



'2k n (+oo,E) 
iv) Proposition 14.11 is proven in the Section [7| 

5 Non— adiabatic Transition Asymptotics 

5.1 The Transition Integral 

From now on, we assume we are in the avoided crossing situation, but we still do not make 
explicit the dependence in the variable 5 > in the notation. 

Section 3 gave us the semiclassical asymptotics of the elements of the S-matrix S(E, e). 
We now compute the small e asymptotics of the integrals that describe the asymptotic states 
i/)j(x, t, e, ±) given by (|4.2jl as \t\ — > oo, for the different channels. 
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We choose our energy density Q(E, e) to be more and more sharply peaked near a specific 
value E G A \ <9A as e —>■ 0. As a result, we obtain semiclassical Born-Oppenheimer states 
that are well localized in phase space. This choice is physically reasonable, and it allows us 
to relate the quantum scattering process to classical quantities. 

More precisely we consider, 



Q(E,e) = e- G(B)/£ e" iJ(B)/£ P(£,e), (5.1) 

where 

CI : The real-valued function G > is in C 3 (A), and has a unique non-degenerate absolute 
minimum value of at E$ in the interior of A. This implies that 

G(E) = g(E-E ) 2 /2 + O(E-E ) 3 , where g > 0. 

C2 : The real- valued function J is in C 3 (A). 

C3 : The complex-valued function P(E,e) is in C 1 (A) and satisfies 



sup 

e>0 



d n 



dE 



< C n , for n = 0, 1. (5.2) 



Remark: Typical interesting choices of Q have G = g (E — E ) 2 , J = 0, and P an 
e-dependent multiple of a smooth function with at most polynomial growth in (E — Eo)/e. 

In our avoided crossing situation, we have already proved the following: A wave packet 
incoming from the left in the remote past produces reflected waves (i.e., components that 
travel to the left in the remote future) that are exponentially smaller than the components 
that travel to the right in the remote future. We have also proved that the non-trivial 
transitions to electronic states that are not involved in the avoided crossing are exponentially 
smaller than those to electronic states that are involved in the avoided crossing. 

Thus, the leading non-adiabatic transitions are described by the asymptotics of those 
coefficients cf (±oo, E, e) that satisfy 

df.(-oo,E, e) = S j<k 5 a - (5.3) 

c-(+oo,E,e) = e~ w ^e l ^ iz ' E)dz/e ' 2 (l + E (e 2 )), (5.4) 

where n = ir(j) = j ± 1. We recall that the error term 0#(e 2 ) depends analytically on the 
energy E in a neighborhood of the compact set A. We have already noted in the comments 
after 1CT23I) that the term E {e 2 ) satisfies (l5~2"j) . 

The form chosen for the energy densities should make it clear that Gaussian wave packets 
will play a particular role in the asymptotic analysis of ()4.5|) . Therefore we use the specific 
notation introduced in (|1.2j) for them. 

Recall that a normalized free Gaussian state propagating in the constant potential e(+oo) 
is characterized by the classical quantities 

A+(t) = A + + itB+, 
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B+(t) = B + , 

a+(t) = a + + r] + t, 

V+(t) = V+, and 

S + {t) = J* (r) 2 + ( S )/2 - e{+oo)) ds, 

with Re (A + B + ) = 1 (see [§]). The associated nuclear wave packet has the form 
e iS +^ 2 MMt), B + (t), e\ a + (t),r] + (t),x) (5.5) 

e «(rr72-e(oo))/ £ 3 j ( _ ( + r]+ t)) 2 B + \ j . 7] + (x - (fl+ + TJ+t)) 1 

ex P "^ o 9 / ^ , -^ n \ r ex P i l 2 [ ' 



7r 1 /4y / e (A + + UB+) ' I 2e 2 (A+ + *££+) J \ e 

We now have everything to state our main result: 

Theorem 5.1 Let ip(x,t,e) be a solution of the Schrodinger equation M.l)) with electronic 
Hamiltonian h(x, S) that satisfies hypotheses H4 ; H5 ; H7. Assume that the solution is 
characterized asymptotically in the past by 

lim \\ip(x,t, e) — ipj(x,t,e, — ) || = 0, 



t— + — OO 

with 



Q(E,e) 



J * J2kj(-oo,E) 



where the energy density is supported on the interval A, and 

Q(E,e)= e~ G ^ 2 e- iJ ^ 2 P(E,e) 
satisfies CI, C2, and C3. Le£ n = IT (J) be given by \3.21 ), and let 



a(E) = G(E) + Im(f k j (z,E)dz), (5.6) 

-'C 

k(E) = J{E) - Re {[ kj{z,E)dz) + u~{oo,E). (5.7) 

Assume E* is the unique absolute minimum of a(-) in Int A. 

Then, there exist So > 0, p > arbitrarily close to 3, and a function eo : (0, So) — > M + , 
such that for all S < So and e < eo{S), the following asymptotics hold as t — > oo . • 

e -i6j(0 £ 3/2 3/4 

1>-(x,t,e,+) = (j) n {x) — 7TU elS+{t)le2 MMt),B + (t),e 2 ,a + (t), V+ (t),x) 

(£a(E(k))\ k " 



x P(E\ e)V¥e- a{Et)/ ^ e -*("(^)-**V(£*))A 2 + ( e -a(e*)AV) + O e (I ft) 
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where tp is parametrized by 



v+ = k* = J2(E*-e n (oo)), a + = k*K'(E*) } B, 



^a(E(k))\ k * 



( d 2 d 2 \ I d 2 

A+ = { ^ a{E{mk ' + l ^< E{ ~ k ^ k * J / y dk^ a{E{k))lk * and (5 - 8) 

S+(t) =t(k* 2 /2-e n (oo)). 

All error terms are estimated in the L 2 (]R) norm, and the estimate 0(e~ a ( E *^ e e p ) is uniform 
in t, whereas O e (l/t) may depend on e. 

Remarks: 

i) All quantities computed from the electronic Hamiltonian h(x, 5) depend on 5, even 
though that dependencs is not specified in the notation. 

ii) The function a has a unique absolute minimum if |A| and 5 are small enough. See 
Proposition 15.11 However, in the case of several absolute minima, one simply adds the con- 
tributions associated with each of them. 

iii) The transitions to states that travel to the left in the future are excluded from our anal- 
ysis because of the lack of uniformity in E in the semiclassical asymptotics of the relevant 
elements of the matrix S(E,e). At the price of some more technicalities, it should also be 
possible to accommodate this situation by our methods. 

iv) When several avoided-crossings are taken into account and meet the requirements of 
Theorem 13.11 c~(oo,E,e) with n = n(j) is given by a product of exponentials of the same 
form as those in ()5.4|) . The analysis of this situation is essentially identical to the single 
avoided crossing situation, mutatis mutandis. 

v) Further properties of ip n (x,t,e,+) are given below. In particular, the characteristics of 
the average momentum k* and its behavior as a function of 6 are detailed in Section 5.2. 
The energy densities corresponding to specific incoming states are studied in Section 6. 
vi) The asymptotics of the incoming wave with the electrons in the state <j)j in the remote 
past are described by the same integral, with the replacements 






2(E-e n (oo)) h+ y/2(E - e^-oo)), 

{ e s ~o. 



Proof of Theorem I5.lt 

Apart from the ^-independent factor given by 

e 



(5.9) 



V2 



*j(0 
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the asymptotics of (|4.5jl are determined by the integral 

P(E,e) 



T(e,x,t) 



a (2(E - e n (oo)))V4 



x e -G(E)/e 2 e -i(tE+J(E))/e 2 e H,(E)/e 2 e i(x^2(E-e n (oc))-u>-(oc,E))/e 2 jfi 

where P(E, e) = P(E, e) (1 + E (e 2 )) satisfies , 
'jj(E) = / kj(z ) E)dz 1 and 

w-(oo,S) = -| o °° (yf2(E - e n (y)) - ^2(E - e n (oo)) ) % 

The (1 + 0e(e 2 )) factor in P(E,e) comes from Theorem 13.11 Recall that 7,,- and u;+(oo, •) 
are analytic in a complex neighborhood of A, and that Im 7j(-E) is a positive, decreasing, 
convex function of E, for 5 sufficiently small. 

In terms of the C 3 functions ()5.6|) and ()5.7|) we can write T(e,x,t) as 

T(e,x,t) = f P{E ' €) 1/A e~ a ^ 2 e -^E + ,(E))^ i^E-etoo))/*? dE 

J a (2(E - e(oo)) 1/4 

where we have dropped the index in the asymptotic eigenvalue e(oo) = e n (oo). In Section [7| 
we analyze the small e asymptotics of T essentially by Laplace's method. The result is 



Lemma 5.1 Let k(E) = \J2{E — e(oo)), or equivalently, E(k) = h e(oo) ; and assume 

that a(-) has a unique absolute minimum E* . For sufficiently small 5, this minimum is 
non-degenerate and satisfies E* G Int A. With k* = k(E*), there exists p > arbitrarily 
close to 3, such that as e — > 0, 

Tie x t) = ^£a(E(k))\ k , +i (x + k*^"(E*)) 
K,,) &a{E{k))\ k .+i(t+& K (E(k))\k-))>' 2 

nE\e) e _ a{E * )/e2 f JtE* + k(E*) - xk*) \ 



x cxp' (x-k-(t + K'{E*)))* Oie-*™*^ 

2e*(£,a(E(k))\ k *+z(t + £,K(E(k))\^ ' 



where the error estimate is in the L (M) norm, uniformly in t. 

Remarks: 

i) If there are several absolute minima, one simply adds their contributions to get the 

asymptotics of T. 

ii) If T is associated with the incoming wave as t — > — oo, the formula holds with Eq in 
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place of E*, k = J2(E — e(— oo)) in place of k*, and the changes 1)5. 9 J) . 

iii) If P satisfies C3 and P(E*, e) = 0(e d ) for some d > |, then the above analysis yields 

no information. 

To relate the integral T to standard Born-Oppenheimer states involving normalized free 
Gaussian states, we must identify (|5.1()jl with ()5.5|) . making use of (|5.8jl . and taking care of 
the x and £ dependence in the non-Gaussian part of (|5.1()|) . That is the content of the next 
lemma which completes the proof of Theorem 15.11 

With the identifications (J5.8J1 . we have 
Lemma 5.2 For small e and < p < 3, we have 

T{e,x,t) = e 3 / 2 2^3/4 ^ (E(fc)) | fc V 1/4 ^^ e -^)A 2 






^q(£(fc))| fc .+i(s + fc*V 
^a(E(fc))| fc . +i(t + ^K(E(k))\ k 

x e^^Vo(A + (£), J B + (t),6 2 ,a + (t),r /+ (t),x) + 0( e -<**>/«V), 
where the error is estimated in the L 2 (M) norm, uniformly in t. 

Furthermore, in the L? norm, for small e and large \t\, we have 

T{e,x,t) = e 3 / 2 2 1 / 2 tt 3 / 4 P(E*,e) Vk* e'< E *^ e -*(«(**)-**V(i5*))A a 

/ -2 \ -i/4 

x e i5 +«/ e2 ^ (A + (t), J B + (t),e 2 ,a + (i),7 7+ (t) ) x) f ^-a(J0(AO)|*.J 

+ 0(e- Q ( B *)/ e2 e p ) + 0( e 3/2 /|£|), 
where the first error term is uniform in t. 

Remarks: 

i) We note that the quantities «(•), fc*, and B + depend only on the index j, while «(•), and 
hence, v4 + depend on both j and n. 

ii) More detailed computations are carried out in the next section, which is devoted to 
specific incoming states. 



5.2 Energy and Momentum Shifts 

When there is a single avoided crossing, we can be more precise about the energy and 
momentum shifts revealed by our general analysis. 

For the rest of this section, we assume h(x, 5) satisfies Hypothesis H6. 
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Under this hypothesis, it is known ^3] that the decay rate in the Landau-Zener formula 
(ET23JI has the form 

T (6) = 6" I (--4) + °^) = 6 " D + <W 3 )> 

and that 

Im^(0,5) = 0(5). 

We use these formulas to get more information on E*, the typical energy of the outgoing 
wave packet, that is determined by the relation 

a'{E*) = G\E*) + Im 7 j(£*) = 0, (5.11) 

where the primes denote derivatives with respect to E. 

In the next proposition, we consider two cases: 

In the first case, we choose the exponent G(E) in the energy density to be independent 
of S. This yields less interesting momentum and energy shifts since they vanish to leading 
order in 5 as 5 — > 0, in keeping with [TU] . 

In the second case, we choose G(E) to depend on 5 in such a way that the incoming wave 
packet contains a sufficiently wide spectrum of energies as 5 — » 0. This implies non-trivial 
behavior of the relevant quantities to leading order in 5. For obvious reasons, we restore 5 
in the notation of this discussion. 



Proposition 5.1 Let 

G(E) = g(E - E ) 2 /2 + 0(E - E )\ 

im ^ E ^ = W) +o(s: ' ) = -m +o(n ani 

a(E,8) = G(E) + Im^(E,S), 

as above. 

i) Assume G is independent of 5. Then, for E* defined by \5.11\) . we have 



E*{8) = E + — ^- + 0(5 



as 5 — > 0. In this case, 



a(E*(5)) = -^P- + 0(5 3 ) > a(E ), and 

kc{E ) 

a"(E*(6)) = g + 0(5 2 ). 
IfG(E) = g(E-E ) 2 /2 + 9l (E - E ) 3 /6 + 0(E - E )\ then 

a{E{6)) = g + 9l ikJ(Eo) +3 kf(Eo) + ° ( 6) - 
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ii) Assume G(E, 5) = L(5(E — Eq)), for some function L, such that 

G(E,5) = g 5 2 (E-E ) 2 /2 + O(5 3 ), 
for some go > 0, uniformly for E e A. Then 

E*(S) = E 1 + 0(5), 
where < E\ — E\{D j g§) is the unique solution to the equation 

(E 1 -E ) = D/(g k 3 c (E 1 )), 
and is independent of 5. In this case, 

u(E*(5)) = 5 2 (-^ + g (E 1 -E ) 2 /2^j + 0(5") 

= S 2 (D^gl'^E, - E Q f" + g (E 1 - E,f/2) + 0(5') 
> a(E ), and 

a"(E*(5)) = 5 2 g + 3^ )+ O(5*). 

Proof: Both statements follow from application of the Implicit Function Theorem and the 

observation that a is a strictly convex function of E on A. I 

Remarks: 

a) The first result shows no effect to leading order in 5 in the exponential decay rate of 
transition probability. The value of E*(5) and the width of the outgoing wave packet can 
be computed. Their variations with respect to the corresponding quantities in the incoming 
wave packet are of order 5 2 , and hence, are rather small. 

b) In case ii) of the proposition, the equation that determines E* can be rewritten as the 
quintic equation 

kl(E) - kl(E)k 2 (E ) - 2D/g = 0. 

c) In the case ii), the variation of exponential decay rate in the transition probability is 
given by 

= a „W W | (^^(M^-M^ + BW)) + o(0- 

d) The results above hold provided one knows E* is the unique absolute minimum of G in 
the set A, which is generically true. Again, if there are several minima, one simply adds the 
corresponding contributions. 
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6 Energy Densities and Transitions when the Incoming 
State has the form (p m 

In this section we study the special case in which the incoming state is asymptotically in the 
past on electronic level j with the nuclear wave function given by one of the functions (p m . 
For simplicity, we restrict attention to wave packets that are incoming from the left. 
In the simplest situation, the incoming wave packet is asmyptotic to 

e ifoi/2-«y(-oo))tA ^ A _ + ^B_,.B_,e 2 ,a_ +^t,r/_,x) <j> s (x), (6.1) 

as t — > — oo. Here we choose ??_ > and the set A, so that rjL/2 + &j{— oo) is in the interior 
of A, and that the minimum of A lies strictly above the spectrum of h(x) for all x. 

We choose a smooth cut-off function F(E) whose support is a subset of the interior of 
A, which takes the value 1 on an interval whose interior contains rjt/2 + ej(—oo), and whose 
length is almost as large as that of A. 

From our assumptions on A, there is a one-to-one correspondence between E G A and 
positive k, such that k 2 /2 + e 3 -(— oo) = E. For E G A, we make the change of variables from 
k to E at t = in the (rescaled) Fourier transform of the Gaussian in 1)6.1)) (see [9 ). Taking 
into account the normalization (J3.7)) of the generalized eigenvectors, this leads to the energy 
density 

Q(E,e) = { } Vo(B-,A_,e 2 , V „,-a_,k(E)) (6.2) 



eJ*J2(E-e j (.-oo)) 



if B_, A_, e% n-, -a-, J2(E - e^(-oo)) 



that we use in (J4.2)) . 

Since r] 2 _/2 + ej(— oo) is in the interior of the set where F(E) = 1, the wave functions (J6.1)) 
and ip(x,t,e, — ) defined by (J4.2J) with the energy density defined by (J6.2J) differ in L 2 (M) 
norm by an 0(e~ c ^ e ) error. To be sure that this error is smaller than the non-adiabatic 
effect we are studying, we assume any one of the following conditions: 

1. Take the avoided crossing gap 5 to be small enough that the non-adiabatic effect is 
larger than the error we make here. 

2. Choose \B_\ to be sufficiently small. That increases the value of C in this error estimate. 

3. Fix the minimum of A, but then choose 7/_ large enough so that the cut off is farther 
out in the tail of the Gaussian in momentum space. This also makes the non-adiabatic effect 
larger since ?]_ is larger. 

With Q(E, e) chosen by (J6.2J) . we have in the notation of 1)5.1)1 . 

U/2(E- ej (-oo)) - ??_) 2 
G{E) = (Re (A./BJ)) -t— J -- — — (6.3) 



\B\- 2 



/D „ (<j2(E-ej(-oo))- 


-v-) 2 


)] 2 




^2(E- ej (-oo))- V _r 





2 
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U2{E - eA-oo)) - r).f 
J(E) = (Im(A./B_)) ^- \ (6.4) 

+ a_ (y/2(E - ej(-oo)) - ri-) - uj(E,-oo), 

P(E,e) = IT- 3 ' 4 e~ 3/2 BI 1/2 (2(E- ej (-oo)))- 1/4 F(E). (6.5) 

Also, conditions CI, and C2 are satisfied, and provided we remove the trivial normalization 
factor of e _3//2 from P(E, e), then condition C3 is also satisfied. 

We already know that asymptotically in the past, the interacting wave function deter- 
mined by ()6.2|) agrees with (j6.1|) up to an 0(e~ c / 6 ) error, and we observe that the density 
Q(E, e) is sharply peaked around the energy 

7/2 



E = — — h ej(-oo) corresponding to rj- = J2(E — ej(-oo)). 



Thus from (|6.3|) . we see that 



G ( E ) = } E ^°L + OUE-Enf). /.<-. <j 



2(V-\B-\y vv u/ n ' •' » (V-\B-\r 

We are not particularly interested in the main component of the wave function for large 
time that has not made a non-adiabatic transition. However, by a similar analysis, it could 
be determined by our techniques. Of course, it is what one would expect from the standard 
time-dependent Born-Oppenheimer approximation. 

Our focus is on the dominant non-adiabatic component, which is determined to leading 
order in e by Theorem 15.11 From the above calculations and Theorem 15.11 we immediately 
get our main result for Gaussian incoming states: 

Theorem 6.1 Assume Hypotheses H4 ; H5, and H7 ; and assume A, A_ ; i?_, a_ ; rj^, 5, 
and the levels j and n have been chosen to satisfy the requirements above. Let ty(x,e,t) be 
the solution to the Schrodinger equation that is asymptotic as t — > — oo to 



^ 2 _/2- ej (-cx)))t/ e 2 



<fo(A- + itB-, 5_, e 2 , o_ + r]-t, ry_, x) 4>j( x )- 



The leading non-adiabatic component of \l/(x, e, t) as t — > oo and e — > in L 2 norm is on 
electronic level <f> n (x) and is given by 

An(e) e W 2 -M°°))*A 2 (po (A + + ttB + ,B + ,e 2 ,a + + r ]+ t,r ]+ ,x) <f> n (x), 

where the values of A +! B + , a + , rj + = k* are those given by \5.§) as in Theorem \5.1\ The 
amplitude for making this transition from level j to level n is given by 



^.( e ) = e ~^ (c) e~ aiE * )/e2 J— e -«E*)-k*a+)/e 2 
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In particular 

-1/2 

■ 4- Tm.^"(E*\k* 2 \ 
v \Bj 2 k* 

l-B-l 



B + = ((G"(E*) + Im^(E*))k* 2 )-^ = (rj^- + Im^(E*)k* 2 



^ + Imi'(E*)\B_\k*' 



(6.6) 



Remark: Depending on the relative size of \B_ | with respect to 5, we can apply Proposition 
15. II to further characterize A n j- 

We now turn our attention to the situation where the incoming nuclear wave packet is 
in the state <p m . The only change from the situation just considered is that we must replace 
the function P(E, e) in ()6.5|) by 

P(E, e) = {-i) m 2~ m ' 2 (ml)- 1 / 2 vr^ 4 e" 3 / 2 (2(E - e j (-oo)))- 1 ' A BI lm+1)/2 (B~) m/2 



x ^( ^ g -y-"- | f(g) . ,, 7) 

Again, this satisfies Condition C3 if we take out the trivial factor of e~^ m+3 / 2 \ 

Theorem 6.2 Assume the Hypotheses of Theorem \6.1\ Let *$?(x, e, t) be the solution to the 
Schrodinger equation that is asymptotic as t — > — oo to 

e ^/2-e j( -oo))t/e 2 ^ m {A_ + UB-, B_, € 2 , a_ + TJ-t, 7]-, x) <j)j{x). 

The leading non-adiabatic component of \I/(x,e, t) as t — > oo and e — > in L 2 norm is on 
electronic level (f> n {x), and is given by 

A%\e) j(vl/2-en(oc))t/e* ^ A+ + itB+j B+ , e 2 , a + + V+ t, V+ , x) <j> n (x), 

where the values of A + , B + , a + , 7] + = k* are those given by h5.8ty . as in Theorem \5.1\ The 
amplitude for making the transition from level j to level n is given by 



m/2 



n > [) " V^- 2-/2( m! )i/2 [ B J nm \e\B.\ 



In particular, B + is again given by 116.6]) and the pre- exponential factor is of order e~ m . 
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7 Technicalities 

Proof of lemma I3.lt We consider only the limit x — > oo and the choice a = +. The 
other cases are similar. In this proof, c n denotes a finite constant that depends only on n, 
but may vary from line to line. 
Explicitly, for any jiGN, 

gn 



y/2(E- ej (x)) = c n (2(E-e J (x))) 1 / 2 - n , 



dE n 
uniformly for E 6 A. So, the first assertion is true. Moreover, 

J^ (y/2(E - ej (x)) - y/2(E - ej (oo))) = 

c n ((2(E - ejix))) 1 ^ - (2(E - e^oo))) 1 / 2 "" 

For n = 0, we have by ()2.2|) . 

2(e i (oo)-e i (x)) 



2(£-e j (x)) - J2(E-ej(oo)) 



2(E-e j (x)) + ^/2(E-e j (oo)) 
0(e i (oo)-e i (x)) = 0(< a; >-< 2+ ">). (7.1) 



For n > 0, we can write 

(2(£ - e^x))) 1 / 2 -" - (2(£ - e^oo))) 1 / 2 "" 

(2(E - e^x))) 1 '* - (2(E - e^oo))) 1 '* 
(2(£-e,-(s)))« 

+ (g (2( g -^g"(c)))^ )^W-^)- ^ 

to which the estimate (|7.1jl applies. The second assertion follows. 
By definition, 

r+(+,x,e) = - £° ((2(E - ej (y))) 1/2 ~ ME ~ ejioo))) 1 /*) dy } 

so that (|7.2jl implies the estimates on Jj^r~j~(+,x,E). 

We now study the properties of the cj's. Again, we shall consider x — > +oo; the other 
case is similar. We first compute 

a^(x,E) = -\ \ = ((<j )j (x), ( f> > l (x))(k j (x,E)+Tok l (x,E)) 

2 Jk j {x,E)k l {x,E)\ 



+ (^-|) {UX \ MX)) ^kfaE) 
{(j) j (x),(l)' l (x)){k j (x,E)+Tak l (x,E)) 



2 ^x~EMx~E) 



+ (-|)%i^*))' (7 ' 3) 
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The presence of the factors (<f)j(x), 4>[(x)) and e[(x), which are independent of E and decay 
as 1/ < x > 2+u , implies together with (J3.17J1 that 

' )il 

C aJ{x,E) = 0(<x>- {2+v) ). (7.4) 



dE n "*J 
We denote the coefficients cj collectively by 



c(x,E,e) = ( C+ S X ^ e l) E <D 



c (x,E,e) 



~i2m 



j 



and the generator of equation (|3.12|) by the 2m x 2m block matrix 
M(x,E,e) = 

e i r o ik ]{ y,E)-k l{ y,E)) d y,^ ^+ ^ R) j f o{k]{ y,E) +h{ y,E)) d y/^ ^-^ R) 

e i f o *(-k ]{ y,E)-k l{ y,E)) d y/S Q -+ { ^ R) j J^-k j{y ,E) +H y,E)) d y^ Q --^ £) j ' 

so that (J3.12J) can be rewritten as 

d 

— c(x,E,e) = M(x,E,e)c(x,E,e). 
ox 

Expressing the solutions as Dyson series, we obtain 

c(x,E,e) = W / •••/ (7-5) 

n=0 Jo Jo Jo 

x M(xi, E, e)M(x 2 , E, e) • • • M(x n , E, t)dx\dx 2 ■ ■ ■ dx n c(0, E, e), 

where, because of ()7.4|) . J °° ||M(y, E, e)\\dy < oo, uniformly for £ 6 A, we get the usual 
bound 

|| c{x, E, e) || < e /o°° IIM(w,s,e)ll* || c (0,S,e)||. 

By showing that \\c(x,E,e) — c(y,E,e)\\ is arbitrarily small for large x and y, we see that 
lim. E _ i . 00 c(x, E, e) = c(oo, E, e) exists. Because of the presence of the phases 
e i/;(rfc J ( 2/) B)+ c rfc ! (3/,s))d 2 //6 2 in M(x,E,e), whose derivatives with respect to E satisfy 

D " jfoirkjiy^+vktiy^dy/e* _ Q(< X > n ) 



dE 
we get for n — 0, 1, 






Thus, by the Lebesgue Dominated Convergence Theorem, we get from (|7.5|) that, as x — > oo 
and uniformly for E G A, 

d n 

c(x,E,e) = 0(1), for n = 0, 1. 
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Finally, we consider 

c(x,E,e) — c(oo,E,e) 



M(y,E,e)c(y,E,e)dy. 



In this expression, we use the above properties of M, c, and their derivatives to obtain the 
last two statements of the lemma. I 



Proof of Lemma 14.11 

We assume t ^ and rewrite the exponential factors in (J4.1J) as 



-t(/ " kj(y,E)dy + tE)/e 



-d-p-^Jo ' k°(y,E)dy + tE)/e 



2 dE 



(* + So ^k](y,E)dy 
Then, for each integral in (|4.1jl . we have 

r Q(E,e)c°(x,E,e) ^ k j iy , E) dy + tE)/S dE 



(7.6) 



kj(x, E) 



xe e 



2 p-^Io k?(y,E)dy + tE)/e 



Q(E,e)c?(x,E,e) 



kj(x,E) (t + f* ^k°(y,E)dy 



E 2 



Ei 







Q(E,e)c?(x,E,e) 



a dE Jkj{x,E) (t + J X ^kJ(y,E)dy 



I € 



2 e ~ i( fo ^(y,E)dy + tE)/e 2 ^ 



(7.7) 



The quantities Cj(x,E,e), Jkj(x,E) } and their derivatives with respect to E are uniformly 
bounded in x and E. Also, 



* dkJ(y,E) 
o dE 



dy 



dk°(±oc,E) 
v — - — 77=. + 0{1) 



dE 



a x 



kj(±oo : E) 



+ 0(1), 



uniformly for £ e A as i ^ icx). 

Thus, the boundary terms in (J7.7J1 satisfy 



xe e 



O 



2 ~i(JokJ(y,E)dy+tE)e 2 



Q(E,e)cUx,E,e) 



kj(x,E) (t + f x ^kJ(y,E)dy 

1 



E 2 



h?(±oo, EJt + x + 0(1] 



+ 



k?(±oo,E 2 )t + x + 0(l) 



(71 



38 



We now apply the restrictions on x/t in the statement of the Lemma. For any choice of 
j and a, they ensure that the denominators on the right hand side of (J7.8J1 can be estimated, 
uniformly in E and for large \x\, by 

JfeJ(±oo, E)t + x + 0(1) | = |x| 1 1 + fcJ(±oo, E)t/x + 0(1/ x 

> \x\(a + 0(1 /\x\)), 

where a is the number that appears in the statement of the lemma. From this, we see that 
the boundary terms in (J7.7J) are 0(l/|x|). 

We estimate the integral term in (|7.7jl in a similar way. Under the restrictions on x/t in 
the lemma, we obtain 



L 



2 e -i($* k°{y,E)dy+tE)/e 



ie e 



^(Q(E,e)(k°(x,E))-V 2 c°(x,E,e)) 
t + Jo ^kJ(y,E)dy 



Q(E, e) (fcjfo E))^ cjjx, E, e) g ^ kj(y, E) dy 
+ Jo ^kJ(y,E)dy 



dE 



= 0(l/\x\). 

This implies the lemma for i ^ 0. When t = 0, the estimate ()7.8|) with £ = yields the 
result in a more direct way. I 



Proof of Proposition 14. H We can write 

ip(x,t,e) - ip(x,t,e,±) = Yl <M X ) x 



(7.9) 



f /■ Q(ig, e) dE e _ itE/e * ^ _ E ^ e _ if: kJ{y , E)dy/e 

{ ^ ^(±00, E) 



+ 



Q(E,e) dE itE/e 2 



c? ±oo,£,e 



a J2kj(±oo,E) 



x e 



-if*k°(y,E)dy/< 



2 _ p -J(xfeJ(±oo,£;)+£j^(±oo,B))/e 2 



+ 



r Q(E,e)c°(x,E,e)dE ^_ UE/e2 kj(±oo, E) - kjjx, E) ^ g fef {y , E )dy/e 2 \ 
J A ylkj{±oo,E)kj(x,E) ^2kj(±oc,E) + ^2kj(x,E) J 



The first step of the proof consists of integrating by parts to get a factor of 1/t according to 



I f(x,E,e)e~ UE ^dE = ^f( x ,E,e)e 

J A t 



itE/e 2 



E-> 







TLmWx-*"** 3 - 



(7.10) 
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We then bound the L 2 (M X ) norm of each term that arises from these integrations by parts, 
with bounds that are uniform in t. 

From the estimates in Lemma 13. 1| we see that all the boundary terms in (J7.10J) coming 
from (|7.9jl are of order < x >~( 1+u > . Thus, their L 2 norms are bounded, uniformly in t. The 
integral terms in (|7.1()j) coming from (J7.9J) all have the form 

/ gJx, E, e) e -Kj^j(y,E) dy +tE)/^ dE ^ ■ = 1? 2^ 3^ (7n) 

J A 

where the first integral from (J7.9J1 contains the function 

m(*,E,e) = {^ ( /=== (^,fi, e) -«?(±°c, £..))) } (7-12) 






(cj(x, £, e) - cJ(±oo, E, e)) £ g^k<J(y, E) dy. 



'kj(±oo,E) 
With the notation of Lemma 13.11 the second integral in (|7.9|) contains the function 



j ( 5( £J ' e ) „<r^„ p ^ ^ ^(i^fiWe^ r # 



cJ(±oo, £, e) (l - e ^(±,^))A 2 ) jf" ^ E) dy _ 



/ k j (±oo,E) 3 V > h dE 

The third integral contains 



kj{±oo,E)kj(x,E) J2k j (±oo,E) + J2k j {x,E) -^ <9£ 



By Lemma \'S. 11 and the condition v > 1/2, each of these functions gj(x, E, e) satisfies the 
following bound, uniformly in E, 

gj(x,E,e) = 0(<x>~ v ) e L 2 (IR). 

Therefore, we can estimate the L 2 norm of the corresponding expression (J7.11J1 by 



J J 9j (x,E,e) e ~i(J * WvWy+tEye* dE 



dx < CAe) 



IM 

where Ci(e) is a finite constant that is independent of t. This finishes the proof. I 

Proof of Proposition 13.11 Since the argument is virtually identical to the one presented 
in ^BJ and ^H], we will be rather sketchy and mainly point out the effects of the parameter 
S and of the non self-adjointness of the generator F(x, S). 
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Expressing the projector P(x, 5) as a integral of the resolvent (F(x, 5) — z) l along a loop 
L (or a finite number of such loops) around the set <j\ (x, S) by means of the Riesz formula, 

P{x,S) = -^- f L (F(x,6)-z)- 1 dz, (7.15) 

we get a bound, uniform in 5 > and x G p a , 

\\P(x,6)\\ < c. 



Indeed, for each x G p a , we can choose the path L uniformly in S by hypothesis. The existence 
of the limits F(±oo, 5) allows us actually to consider only a finite number of distinct loops 
a finite distance g/A away from spectrum of F(x, 5), for all (x, 5) . Also, uniformly in 5 > 0, 

\\(F(x,S)-z)- 1 \\ < c, (7.16) 

for z on the corresponding loop L, since | det(F(x, 5) — z)\ > ((?/4) n and F(x, 5) is uniformly 
bounded. By a similar argument, using hypothesis H4, we get, uniformly in 5 

\\P(x,8)-P(±oo,S)\\ < 



as x — > ±oo in p a . With the notation ' for Jr-, we get from ()7.15jl . 



< x > 2+u ' 

i,ro rrof 



P'(x,S) = — i (F(x,6) - z)- 1 F'{x,S) (F(x,6) - z)' 1 dz. 
Thus, hypothesis H4 yields, uniformly in 5, 

\\P'(x,8)\\ < c -—- % 

and a similar uniform estimate for K(x, 5) = [P'(x, 5), P(x, 5)], 

\\K(x,6)\\ < T ^,. (7.17) 

The operator K is the generator of the intertwining operator W defined by 

W'(x, x ,S) = K(x,S)W(x,x ,S), with W(x ,x ,5) = E 
It satisfies 

W(x,x ,5)P(x ,5) = P(x,S)W(x,x ,S), (7.18) 

for all (a;, 8) (including x = ±oo). 

Following [TU] . we construct a hierarchy of generators. Let F (x, 5) = F(x, 5), 
Po(x, 5) = P(x, 6), and K (x, 5) = K(x, 8). For q G N*, we inductively define 

F q (x,S,e) = F(x,5) - e 2 K^^S^e), 
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assuming e is small enough so that the spectrum of F q is separated into two disjoint parts 
corresponding to those of F. We define P q (x, 5, e) to be the spectral projector for F q corre- 
sponding to P(x, 5) as e — > by perturbation theory. Then, 

K q (x,8,e) = [P q (x,8,e),P q (x,8,e)}. 

Sections II. A and II. B of (THj and ()7.16|) and (J7.17J1 show that there exist constants e* > 0, 
r > 0, T > and c > 0, all independent of 5 > 0, such that for all e < e*, all x G M, and 

q = q* = [r/e 2 ], 

||^*-iM,e)|| < <x ° >2+ ^ (7-19) 

e -r/, 2 
\\K q *(x,5,6)-K q *- 1 (x,5,e)\\ < c - - >2+i/ . 

We define 

F*{x,8,e) = F q *{x,8,e) = F{x,8)-eK^ x {x,8,e), (7-20) 

P*(x,5,e) = P q *(x,8,e), 
K*(x,6,e) = K q *(x,S,e), 

and the evolution operators W* and S* by 

Wl(x, Xo, 8, e) = K*(x,5,e) W*(x, Xo,S, e), with W(x 0l x ,5, e) = iZJ 

and 2e 2 S' + (x,x ,(5, e) = W*(x ,a;, #> e ) ^*(^j £, e) W*(x, x , 6, e) S*(x, a? 0) 8, e), 

with S*(x ,a;o,(5, e) = if (7-21) 

The intertwining property (|7.18jl still holds with the * indices. Therefore, S* satisfies 

[E*(x,Xo : 8,e), P*(xo,8,e)] = 0, for all x E 1R. 

It follows from the definitions that the operator 

V*(x,x ,8,e) = W*(x,x Q ,8,e)E*(x,x Q ,6,e) 

satisfies 

ie 2 V*(x,x ,d,e) = {F*(x, 5, e) + ie 2 K*(x, 5, e)) V*(x, x , 5, e) 

and 

V*(x,xo,6, e) P*(x ,5, e) = P*(x, 5, e) K(x,^o,^,e). (7.22) 

Moreover, 

^7 e (x, x ,5) - V*(x,x ,8,e) = (7.23) 

* / V*(xo,y,5,e)(Kf(y,5,e)-K <r -i(y,5,e))U e (y,xo,8)dy. 
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The proposition will follow from 

U e (x,x ,S) - V»(x,x ,6,e) = 0(e~ v/e2 ) } (7.24) 

$Hty and 

lim PJx, S, e) = P(±oo,6), 

x — >±oo 

due to (I7~TT31) and (JH|. 

To prove (|7.24|) . we first prove that V* is uniformly bounded in x, Xo, and e. The analysis 
leading to Lemma 13.11 implies that U e is uniformly bounded in x, Xo, and e. This property 
is a consequence on the fact that the eigenvalues of F are simple and real, so that the 
decomposition (J3.11J) holds and the singular exponential factors are phases. Note that the 
lack of orthogonality of the eigenprojectors of F(x, 5) makes the bound on U e dependent on 
8. 

Choose B(5) > 0, such that ||?7 e (x, Xo, S)\\ < B(S). From (J7.23J1 we get the inequality 

\\V*(x,x ,5,e)\\ < B(8) 1 + / sup \\V*(y,y ,5,e)\\ dx, 

\ Jxo 2/0,2/ \ x I ) 

for some C. This implies that for some C, the quantity v(e,S) = sup ||K(aj, xo, S, e)|| 

x ,x 

satisfies 

v{e,8) < B(8) (l + Cv(e,5)e- r/e2 ). 

This implies 

v(e,S) < 3S < v(S), 

where v(S) is uniformly bounded for sufficiently small e. 
We now use (J7.23J) again to see that 



U e (x,x ,S) - V*(x,x ,6,e) || < v{8) / i „\ 2+v — rfx 



g((5)C e- r / £ 



This proves ()7.24|) and completes the proof of the proposition. I 

Proof of Lemma I3.3t Degenerate perturbation theory for self-adjoint matrices and 
hypothesis H5 (see [13]) show that there exist f(z,5) and p(z,8), analytic in z for fixed S, 
and C l as functions of (z, 8), such that 



ej (z,5) = f(z,5) - \jfa$) ( 7 - 25 ) 

e n (z,S) = f{z,5) + - ^Jp{z,5). 
where, as (z,5) — > (0,0), 
f(z,6) = f(0,0) + O(\z\+S) = e c + 0(\z\+5), and p(z,S) = 0(\z\ 2 + 5 2 ). 
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Moreover, p(z,8) has two simple zeros, the complex crossing points, Zq(8) and z Q (8) that 
have Zq (8) = 0(5). For concreteness, we arbitrarily choose e 3 - < e n on the real axis, although 
this is irrelevant for the analysis. Thus, by H7, we can write 



p(E-e 3 (z,8)) = y/2(E-f(z,6)) 1 + 



>Jp(z, <*) 



1/2 



2(E-f(z,8)) 
where (E — f(z, 5)) and its inverse are analytic in p a , uniformly in E G A. Moreover, 



Therefore, since J 2 (E — f(z,8)) is analytic, and we can choose the loop ( encircling zq(8) 
or z (8) to satisfy \(\ = 0(8), we see that 

f ^2(E-e 3 (z,8))dz = I J X^ dz + 0(8"), 

J ( ^2(E-e c ) J C * 

and 

^<fe = 0(f). 

C 2 

r Jp(z,8) r 

In these two expressions, / - — dz — e 3 (z, 8) dz due to the analyticity of / in 

(J7.25J) . Taking the imaginary part yields the first statement of the lemma. Note that 
we do not have to worry about sign issues because Theorem 13.11 ensures the decay rate, 



Im / \J2(E — e 3 (z, 8)) dz, is positive. 



The two other statements follow from similar considerations for the integrals 

— — Im / \ 2{E — eAz,8)) dz = Im / . dz and 

dE k V l A " k p(E-e 3 (z,8)) 

d 2 r I r 1 

-—-rim / \2(E - ej(z,8)) dz = - Im / ^ dz. 

dE 2 i C V V J k (2(E-e 3 (z,8))f 2 



Proof of Lemma 15.11 

Consider first the minimization of the negative of the real part of the exponent. 

Since 'jj(E) tends to zero with 8 (absent in the notation), if 8 is small enough, we must 
look for minima in a neighborhood of Eq that satisfy the equation 

a'(E) = g(E-E ) + Im ^(E) + 0(E - E ) 2 = 0. 
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We consider the absolute minimum E* of a and assume it is unique. By Lemma 
Im jj(E) < 0, so E* > E Q . Note that E* does not depend on x or t. Also, Im ^(E) = 0(8), 
n = 0, 1, 2, uniformly in E. So, we can assume E* is non-degenerate since 

a"(E*) = g + Im i-{E*) + 0(E* - E ) > 0. 

In terms of the variable k e [fci, fc 2 ], we view T as the (scaled) inverse Fourier transform 
of the function 

R(k, t, e) = V^ 2 e - a{E{k))/e2 P(E(k), e) Vk e -Mm)l* e -^ 2 /2+e(oo))/^ ^^(/^ 
where xs(") is the characteristic function of the set S. That is 

T(x,t,e) = (JT 1 i?(^,e))(x), 
where JF £ is defined by 

(F e g)(x) = -=L= / g(k)e-^ 2 dk. 

With the variable k G [fei, fa] we have 

<9 2 



dk 2 
and expanding around /c*, 
T(e,x,t) = e"^*^ 2 



a(E(k)) 



k* z a"(E*), 



x f Vk P(E(k), e) e- ~ g * (fc - fc,)2 e « fe - fc *) 3 )A 2 ^(M,*)/^ dife, 

•Afcl,fc 2 ] 



where the negative of the imaginary part of the exponent is denoted by 

P(k,x,t) = t ( he 2 (oo)j + K(E(k)) - xk. 

We now introduce fi(e) = e s > 0, with 2/3 < s < 1. It goes to zero in such a way that 

H{e)/e » 1 and /i(e) 3 /e 2 << 1. 

Because E* is a unique absolute minimum, the behavior of at(E) close to E*, and the 
assumption (|5.2|) on P, we can reduce the integration range in T to [fc* — /i(e), fc* + A*( 6 )] a ^ 
the expense of a relative error whose L 2 norm is of order 0(e°°), uniformly t. More precisely, 

T(x,t,e) = (^- 1 (R 1 + R 2 ))(-,t,e)Kx), 

where 

Ri(k,t,e) = X\k*-^e),k*+fi(e)){k) R{k,t,e), and 

R 2 (k,t,e) = X[k*-n{e),k*+n{e)]c{k) R{k,t,e). 
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For some a* > and r > 0, 

\R 2 {k,t,e)\ < r e - a{E * )/e2 e~ a * Me)/e)2 e \VkP(E(k),e)\. 
Hence, by the Parseval identity, uniformly in t, we have 

||^T 1 (i2 a )(-,t,e))|| = (/ \R 2 (k,t,e)\ 2 dk) = 0( e -^*)A 2 e°°) 



'[**- M (e),fc*+M(e)]° 

In the remaining integral containing Ri, we further estimate 

e o{k-k*fh 2 = l + 0(/i(ef/e 2 ) = l + 0(e 3s - 2 ), (7.26) 

and 

VkP{E{k),e) = Vk? P(E*, e) + 0(ji(e)) = V¥ P{E*,e) + 0{e s + e 2 ). 

The contribution of order e 2 comes from the error in the computation of the coefficient c~. 
Using the Parseval identity again with uniform bounds on the exponential factors of R\, we 
see that the contribution to T coming from the error term 0(e s ) is bounded uniformly in t 
in the L 2 (IR X ) norm by 0(e~ a ^ E *>' e e 1+2s ). Similarly, the error term stemming from ()7.26|) 
yields an error in the L 2 (IR X ) norm of order 0{e~ a ^ E '' e e 4s_1 ). 
To compute the leading term, we expand f3(-,x,t) around k* as 

p(k,x,t) = tE* + k(E*) - xk* 

+ (k-k*) (k*t + ^r<E{k)) - x\ 

h- ( ^f- m^t- (7 - 27) 

where k lies between fc and k*, and the third derivative is independent of t and x. The last 
term in ()7.27|) gives rise to a contribution which is of order 0(e~ a ( E *" e e 4s_1 ) in the L 2 (IR X ) 
norm, uniformly in t, as above. 
Therefore, in the L 2 sense, 

a{E*)/e 2 -i(tE* +n(E*)-xk*) / e 2 



T(e,x,t) = e~ a[E )/e e 



[ Vk*P(E*, e) g-^^'i+I^Ef*))!*.)-^) 

J\k*-u(e),k*+u(e)} 



X 

l{k*-n(e),k*+fi(e)] 

x e - i ^^ aiEmk * +i{t+ -& KiEm ^ ) dk + 0(e p ) + 0{e°°) 

where p = min(l + 2s, 4s — 1) G (0, 3) can be chosen arbitrarily close to 3. Again, at the 
cost of an error whose L 2 norm is 0(e~ a ^ E *^ e e°°), uniformly in t, we can extend the interval 
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of integration to the whole real line and compute the Gaussian integral explicitly according 
to the formula (for Re M > 0) 



k * e -(M(k-k*f /2+iN(k-k*))/e 2 dk 



e 2e-* M 




(k* - iN/M) 



We then get the result with 
d 2 



M 

N 



dk2 <E{k))\ k ,+i{t+^ 2 K{E{k))\ k ^ 

k*t + -?-K(E(k))\ k * - X. 
ok 



and 



Proof of Lemma 15. 21 The first assertion is straightforward. The second follows from 
the identity 

<Po(A+(t), B+(t), e 2 , a+(t),7i+(t), x)x = y {A + {t), B+(t), e 2 , a+(t), r] + (t),x) (x - a+(t)) 

+ v?o04+(£), #+(£), a + (£),e 2 ,r/ + (£), a;,) a+(t). 

The first term is 0(e) in L 2 (M) by scaling, and the second is of order a + (£) = k*t(l+0(l/\t\)) 
for \t\ large. We insert this in the first part of the lemma to obtain the second part as t — > ±oo. 
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